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ABSTRACT 

New optical spectroscopy of the HMXB microquasar LS I +61 303 is presented. Eccen- 
tric orbital fits to our radial velocity measurements yield updated orbital parameters in 
good agreement with previous work. Our orbital solution indicates that the periastron 
passage occurs at radio phase 0.23 and the X-ray/radio outbursts are triggered 2.5-4 
days after the compact star passage. The spectrum of the optical star is consistent 
with a B0 V spectral type and contributes ~65 percent of the total light, the remain- 
ing due to emission by a circumstellar disc. We also measure the projected rotational 
velocity to be usini ~ 113 km s _1 . 
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1 INTRODUCTION 

LS I +61 303 is a peculiar radio and X-ray star which 
has attracted the interest of astronomers during the last 
quarter of a cent u ry. A ttention to it was first called by 
IGregorv fc Tavlorl il9T8l) during a galactic plane survey 
aimed at the identification of radio variables on both short 
and long time-scales. This precursor work unveiled what has 
remained a key property of LS I +61 303 over the years, that 
is, a clearl y periodic modulation o f its radio emission every 
26.5 days (iTavlor fc Gregorvlll982ft . 

Just to mention the key facts, LS I +61 303 is a high 
mass X-ray binary (HMXB) with the optical component be- 
ing a B0 Ve star with a circumstellar disc, V ~ 10.7 mag, and 
with a n estimated distance of 2 kp c (lHjrtchin£s_Ji_Cj^inr>tod 
Il98ll: iParedes fc Figuerasl Il986t iFrail fc Hiellmind ll99llL 
The periodic modulation has also been observed at op- 
tical and infrared w avelengths (lMeridds£n_j£_Mazehl Iroiit 
IParedes et"aflll994h . in X-rays JParedes et alJll997TI . and 
possibly in 7-rays as well jMassill2004air This modulation is 
widely attributed to the orbital period of the binary system 
and the most accurate determinati on (P = 26.49 60 + 0.0028 
d) comes from radio observations iGregorvl |2b02) . A second 
modulati on, on a longer time scale of ^4 yr ; has also been 
reported iParedeslll987l : IGregorv et al.lll989t) and is likely 
associated with an outward-moving density enhancement or 
shell in the equatorial d isc of the rapidly rotating Be star 
iGregorv fc Neisrj l2002V The recent detection of relativis- 
tic radio jets - at milli- arcsecond angular scales and with 
hints of precession - bv lMassi et alJ (I2004T) firmly includes 



LS I +61 303 in the selected class of galactic microquasars. 
The nature of the compact object remains unknown given 
the poor constraints in the system parameters and the lack 
of classical neutron star signatures (e.g. X-ray pulses or Type 
I bursts). 

The seminal paper bv lGregorv fc Tavloil lll978l) already 
pointed out the possible connection of the peculiar radio 
star with the COS B 7-ray source CG135+01. The idea 
of LS I +61 303 as a 7-ray emitter has persisted over 
the years. Today, this star is the best counterpart can- 
didate to 3EG J0241+6103, one of the unidentified high 
energy (> 100 MeV) sources in the 3rd EGRET catalog 
jHartman et alJll999T) . Such association is likely not to be a 
unique case. The microquasar LS 5039, probably associated 
with another EGRET source, brought the idea that micro- 
quasars could be behind som e of the unidentified EGRET 
sources (IParedes et alJ koOOV Theoretical modelling of the 
7-ray emission strongly supports this interpretation. Inverse 
Compton scattering of stellar photons by the relativistic 
electrons in the jets, can account f or the overall properties of 
micro quasars from radio to 7-rays iBosch-Ramon fc Parede"s! 
2004). 



In this context, it is surprising that the orbital param- 
eters of LS I +61 303 have remai ned poorly known over 
the ye ars. Over two decades ago, iHutchings fc Gramotonl 
( 1981) carried out a series of optical spectroscopic observa- 
tions from which a radial velocity curve of the binary system 
was determined. Spectroscopic determination of the radial 
velocity curve of LSI+61 303 is complicated by contamina- 
tion of the narrow photospheric lines by emission from the 
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Be circumstellar envelope. Here we present new optical spec- 
troscopic observations of LS I +61 303 using modern tele- 
scopes and detectors. The idea was to obtain an improved 
set of spectroscopic orbital elements which could be used to 
address issues such as the system's mass function, orbital 
eccentricity, or the phase of periastron, among others. An 
accurate knowledge of such fundamental properties is es- 
ential to remove the lingering uncertainties that plague all 
theoretical models for the multi-wavelength variable emis- 
sion of LS I +61 303, which we believe is a prime target to 
attain a better understanding of the 7-ray sky. 

In the following sections we present our observation and 
reduction procedures, orbital fits and discussion of results. 



2 OBSERVATION AND DATA REDUCTION 

LS I +61 303 was observed using three different telescopes 
and instruments in the blue spectral range (AA3400-5800) 
during years 2002 and 2003. Seventeen intermediate reso- 
lution spectra (0.63 A pix" 1 ) were obtained using the In- 
termediate Dispersion Spectrograph (IDS) attached to the 
2.5 m Isaac Newton Telescope (INT) at the Observatorio 
del Roque de Los Muchachos on the nights of 27-31 July 
2002, 13 March 2003 and 1-7 July 2003. We used the R900V 
grating in combination with the 235 mm camera and a 1.2" 
slit to provide a spectral resolution of 83 km s" 1 (FWHM). 
The spectral type standard r Sco was also observed with the 
same instrumentation for the purpose of radial velocity anal- 
ysis. In addition, this star has a very low p rojected rotational 
veloci ty, low er than 5 km s -1 ac c ording to lHardorp fc Schold 
Jl970h and IStickland fc Llovdl lll995l) . Six high-resolution 
spectra of LS I +61 303 were also obtained with the 2.7 
m Harlan J. Smi th telescope and the 2dcoude spectrograph 
jTull et alJll995h at McDonald Observatory (Texas) on UT 
3-5 August 2002 and 8-9 April 2003. We used the F3 focal 
station, the E2 echelle grating, and a l"2 slit to achieve a 
FWHM resolution of about 5 km s _1 . Finally, five low res- 
olution spectra (2.74 A pix -1 ) were obtained on the nights 
of 13 October and 15-17 December 2003 with the ALFOSC 
spectrograph at the Nordic Optical Telescope (NOT) at the 
Observatorio del Roque de Los Muchachos. We employed 
grism 6 and a 0.5" slit which resulted in a spectral resolu- 
tion of 235 km s _1 , as measured from gaussian fits to the 
arc lines. A complete log of the observations is presented in 
Table [U 

The images were bias corrected and flat-fielded, and 
the spectra subsequently extracted using conventional op- 
timal extraction techniques in order to o ptimize the signal- 
to- noise ratio of the output (Home 1986). Frequent observa- 
tions of comparison arc lamp or hollow cathode lamp images 
were performed in the course of each run and the pixel-to- 
wavelength scale was derived through polynomial fits to a 
large number of identified reference lines. The final rms scat- 
ter of the fit was always <l/30 of the spectral dispersion. 



3 RADIAL VELOCITIES AND ORBITAL 
SOLUTION 

All the spectra were prepared for the cross-correlation anal- 
ysis by subtracting a low order spline fit to the continuum 



Table 1. Log of the observations. 



Date 


Object 


Obs. 


Exp. Time 


Dispersion 








(seconds) 


(A pix- 1 ) 


2002-07-27 


LS I +61 303 


INT 


600,300 


0.63 


2002-U i-zo 


TOT K1 Q n Q 
LD 1 + Ol 303 


INT 


2x< >00 


0.63 


onno n*7 on 
2002-0 i-Z\) 


tot 1 a-\ n 
Ld 1 +61 303 


INT 


nr\r\ ocn 
600,300 


0.63 


nnnn at r i(\ 
2002-0 /-OO 


to t 1 fti n 
LD 1 +61 303 


INT 


2x000 


0.63 


nnnn n*7 01 
2002-0 /-31 


tot 1 a-\ n 
Ld 1 +61 303 


INT 


30(1 


0.63 


2002-08-03 


LS I +61 303 


McD 


600 


0.11 


2002-08-04 


LS I +61 303 


McD 


800 


0.11 


2002-08-05 


LS I +61 303 


McD 


600 


0.11 


2003-03-13 


LS I +61 303 


INT 


2x300 


0.63 


2003-04-08 


LS I +61 303 


McD 


1200 


0.11 


2003-04-09 


LS I +61 303 


McD 


2x1200 


0.11 


2003-07-01 


LS I +61 303 


INT 


300 


0.63 


2003-07-01 


T SCO 


INT 


1 


0.63 


2003-07-02 


LS I +61 303 


INT 


500 


0.63 


2003-07-03 


LS I +61 303 


INT 


500 


0.63 


2003-07-04 


LS I +61 303 


INT 


(iOO 


0.63 


2003-07-05 


LS I +61 303 


INT 


400 


0.63 


2003-07-06 


LS I +61 303 


INT 


500 


0.63 


2003-07-07 


LS I +61 303 


INT 


750 


0.63 


2003-12-15 


LS I +61 303 


NOT 


60 


2.74 


2003-12-16 


LS I +61 303 


NOT 


60 


2.74 


2003-12-17 


LS I +61 303 


NOT 


60 


2.74 



and subsequently rebinned into a uniform velocity scale of 
42 km s . The template spectrum, r Sco, was broadened 
by 113 km s -1 and 226 km s" 1 in order to match the width 
of the photospheric lines in the LS I +61 303 spectra ob- 
tained with the INT and NOT, respectively. These broad- 
ening factors were derived by subtracting several broadened 
versions of r Sco from individual spectra of LS I +61 303 
and performing a y 2 test on t he residuals (see details in 
iMarsh. Robinson fc Woodlll994h . To broaden t he template 
spectrum we employed a Gray rotational profile jGravll992l) 
with a limb darkening coefficient of e = 0.33, appropriate for 
our spectral type and wavelength range. Note that the opti- 
mal broadening is driven by differences in instrumental reso- 
lution for the case of the NOT spectra, but not the INT since 
the same setup was used in both the LS I +61 303 spectra 
and the template. In the latter case, the broadening is di- 
rectly associated to v sin i and confirms that the optical star 
in LS I +61 303 is rapidly rotating. On the other hand, since 
the instrumental resolution of the McDonald data is much 
higher than the template, we decided to optimally broaden 
the target spectra (to match the template) for the sake of 
the cross-correlation analysis. The optimal broadening was 
found to be 113 km s - . 

We cross-correlated every individual spectrum of 
LS I +61 303 with the conveniently broadened template of 
t Sco after masking out the interstellar absorption features 
at A3934 (Ca II K) and A4430. The correlations include H I 
lines from H/3 to Hg and all He I and He n lines in the spec- 
tral range AA3850-5020 (see Fig. 3). Furthermore, since the 
H I lines are prone to contamination by stellar wind emis- 
sion in early-type stars, we also decided to obtain a set of 
radial velocities with only He I and He n lines. T he result- 
ing ve locities were folded on the radio ephemeris of lGregorvl 
(2002) (P = 26.4960 d, T =HJD2443366.775). TableHlists 
the radial velocities measured using the two different sets of 
spectral lines. Individual velocities were extracted following 
the method of Tonry & Davis (1979), where parabolic fits 
where performed to the peak of the cross-correlation func- 
tions, and the uncertainties are purely statistical. 
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Table 2. Radial velocities of LS I 
correlation. 



-61 303 measured from cross- 



Table 3. Orbital solutions. 



HJD- 


Phase 


RV (H+Ho) (km s' 1 ) 


RV (He) (km s _1 ) 


Obs. 


2452000 




Value a 


Value a 





483.7265 


0.0878 


-48.3 


2.3 


-37.6 


2.7 


INT 


483.7336 


0.0881 


—41.2 


3.3 


— 29.5 


4. 1 


INT 


484.7173 


0.1252 


—41.1 


4.5 


—34.6 


5.5 


INT 


484.7250 


0.1255 


—44.3 


2.9 


-37.2 


3.5 


INT 


485.6677 


0.1611 


— 31.0 


2.4 


— 22.1 


2.9 


INT 


485.6735 


0.1613 


— 32.2 


4.9 


— 20.5 


fi.O 


INT 


486.6493 


0.1982 


— 22.5 


2.3 


— 11.4 


2.(i 


INT 


486.6593 


0.1985 


-17.1 


2.1 


-8.2 


2.5 


INT 


487.6722 


0.2368 


-15.6 


3.0 


-4.1 


3.5 


INT 


489.9708 


0.3235 


-54.5 


3.0 


-42.2 


4.2 


McD 


490.9739 


0.3614 


-47.4 


0.8 


-42.3 


1.0 


McD 


491.9774 


0.3992 


-56.2 


1.5 


-64.9 


1.9 


McD 


712.3223 


0.7154 


-57.8 


9.4 


-41.8 


11.6 


INT 


712.3278 


0.7156 


-70.1 


6.3 


-59.0 


8.0 


INT 


737.6144 


0.6700 


-60.0 


1.1 


-40.0 


1.1 


McD 


738.5838 


0.7066 


-49.0 


1.3 


-41.9 


1.(1 


McD 


738.6003 


0.7072 


-40.1 


1.2 


-36.7 


1.4 


McD 


822.7259 


0.8822 


-73.1 


2.0 


-64.9 


2.5 


INT 


823.7330 


0.9202 


-58.2 


1.9 


-49.9 


2.3 


INT 


824.7324 


0.9579 


-56.7 


4.4 


-51.2 


5.4 


INT 


825.7301 


0.9956 


-45.2 


2.4 


-34.4 


2.9 


INT 


826.7322 


0.0334 


-52.2 


2.4 


-43.6 


2.9 


INT 


827.7361 


0.0713 


-36.3 


1.9 


-28.7 


2.3 


INT 


828.7309 


0.1088 


-32.0 


6.9 


-26.0 


8.3 


INT 


989.3004 


0.1690 


-31.5 


6.3 


-12.1 


8.2 


NOT 


990.2950 


0.2065 


-29.6 


4.8 


-21.3 


5.8 


NOT 


991.2943 


0.2442 


-38.1 


4.9 


-23.5 


6.0 


NOT 



The radial velocity curve was modelled with an eccen- 
tric orbital solution. This was accomplished using the so- 
phisticated binary fitting program originally developed by 
IWilson fc Devinnevl Jl97lt hereafter W-D) 1 . Briefly, the W- 
D program assumes a Roche model to describe the shape of 
the components and includes a complete treatment of the 
radiative properties of the system (limb darkening, temper- 
ature variations accross the surface, etc.). Thus, the model 
adequately takes into account the non-keplerian contribu- 
tions to the measured radial velocities. Fitting is carried out 
through differential corrections and, in the case of radial 
velocity solutions, convergence is achieved very rapidly. In 
the solutions we adjusted the following orbital parameters: 
velocity semiamplitude (K\), systemic velocity (7), eccen- 
tricity (e), argument of the periastron (w), and phase of the 
periastron (<^ P eri)- A set of fixed parameters necessary to the 
fitting procedure were tuned to match the expected radiative 
properties o f the B0 V optical companion (i.e., T e ff«28 000 
K, log a~4: ICoxll200oll . although the effect of the adopted 
parameters on the solutions is very weak. Table |3] presents 
the final best-fitting parameters for the two sets of radial 
velocities. 

The two solutions yield orbital parameters that are con- 
sistent within la except for the systemic velocity, which is 
more negative when the Balmer lines are included in the 
correlation. This suggests that the Balmer lines may be con- 
taminated by some blueshifted absorption, consistent with 
the suspicion that the Be primary might have a significant 
wind. Therefore, we prefer to adopt the solution that con- 

1 The original code has suffered major upgrades since its first 
release, including significant improvements in the underlying 
physical models. The most recent version of the W-D pro- 
gram together with the relevant documentation can be found in 
ftp : //astro . uf 1 . edu/pub/wilson/lcdc2003 



Parameter 


H i+He i+He 11 


He I+He II (adopted) 


c 


0.63±0.11 


0.72±0.15 


u (deg) 


7.8+13.9 


21.0±12.7 


7 (km s — 1 ) 


-48.4±1.8 


-40.2±1.9 


'/'peri 


0.22±0.02 


0.23+0.02 


Xi (km s" 1 ) 


18.1±3.3 


22.6±6.3 


ai sini (Rq) 
f(M) (Mq) 
a (km s _1 ) 


7.4±1.6 
0076+ 0060 

u - uu '"-0.0038 

7.8 


8.2+2.9 
0107+ - 0163 

u ' uiu ' -0.0077 

8.5 



\ -; 
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• 
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Figure 1. Radial velocity curve obtained by cross-correlating 
with t Sco using the He I and He II lines in our spectral range. 
The best-fitting solution, using an eccentric orbit, is overplotted. 
The upper panel shows the residuals of the fit. 



siders the velocities derived from the He lines alone, which 
should be largely free from contaminating emission. Figure^ 
shows these velocities, together with the best-fitting orbital 
solution. To better illustrate the resulting orbit, we show 
in Fig. [5] the relative motion of the compact object around 
the optical companion as seen from above (i.e., no inclina- 
tion value assumed). Some relevant phases (i.e., periastron, 
apastron, and conjunctions) are labelled. 

We can now use our orbital solution to Doppler-shift all 
the individual spectra and produce an average spectrum in 
the reference frame of the primary (see Fig. E| . This can be 
compared to the spectrum of the B0 V template r Sco and 
we find that LS I +61 303 is best fitted when the template is 
scaled by a factor 0.65. In other words, the results indicate 
that the primary contributes 65 percent of the total light. 
The rest can possibly be attributed to emission from the 
disc around the Be star. The residual of the spectral sub- 
traction indeed shows broad double-peaked Balmer emission 
features (i.e., characteristic of a disc geometry) with possi- 
ble narrow reversal absorptions superimp osed, in accordance 
with observations of Be-shell stars fe.g.. IPorter fc Riviniusl 
l2003t ISlettebak1ll98ct . However, inspection of high resolu- 
tion Ha profiles published in literature (e.g. Paredes et al. 
1994) show clear double-peaked profiles with shallow central 
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Figure 2. Relative orbit of the compact object around the optical 
component, which lies in the ellipse focus at (0,0). Since the actual 
scale of the orbit is not known (both the inclination and the mass 
ratio cannot be determined with the present single-lined radial 
velocity data), the coordinates are in units of the orbital semi- 
major axis. Relevant phases such as the periastron, apastron, and 
conjunctions arc indicated. Stars mark 0.1-phase intervals. 
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Figure 3. Doppler corrected average of LS I +61 303 in the 
rest frame of the optical star together with the B0 V template r 
Sco and the residual after an optimal subtraction (top) . Features 
labeled "IS" are of insterstellar origin. The template has been 
broadened by 113 km s — 1 and scaled by a factor 0.65 to match 
the depth of the absorption lines in LS I +61 303. 



absorptions, reminiscent of the intermediate inclination Be 
case proposed by Hanuschik (1996). Therefore, contamina- 
tion of the residual Balmer profiles by shell lines can almost 
be ruled out and we interpreted the observed profiles as be- 
ing shaped by kinematic (doppler) broadening. 



4 DISCUSSION 

Our radial velocity data and orbital solution yields a num- 
ber of orbital parameters but does not provide a full charac- 



Figure 4. Mass constraints for the two stars in LS I +61 303 de- 
rived from our orbital solution. The lines have been calculated for 
different inclination values and using our adopted mass function 
(top), plus (bottom left) and minus (bottom right) one sigma. 
The shaded region represents the interval of likely masses for the 
B0 V optical companion (see text). 



terization of the system's physical properties. For example, 
only the mass function and not the mass ratio can be deter- 
mined because only the velocities of the optical component 
are measurable. In addition, the orbital inclination remains 
unknown as usual in any radial velocity curve solution. Thus, 
with the available data we can only plot Fig. 2] that relates 
the mass of the compact object (Mx) and the optical com- 
panion (M opt ) for several values of the orbital inclination. 
Note that the plot has been computed with our adopted ra- 
dial velocity solution that yields /(M)=0.0107t5}:oo77 M©. 
To further constrain the system's properties we must employ 
additional external information. 

The expected masses of the components can be used 
to constrain the system's geometry. For example, the mass 
of the optical component of the system can be inferred 
from its atmospheric properties. Discussions on the most 
li kely mass values for the primary star can be foun d 
inlHutchings fc Cramptonl il98ltt. iMartf fc Paredesl il995t) . 
iPunslvl lll999l) . and lMassi et alTll2004l) . overall favouring a 
mass interval of M op t =10-15 Mq (shaded area in Fig. Hi. 
With this constraint, the compact object would be a neu- 
tron star for inclinations 25° < i < 60° and a black hole if 
i<25°. 

There are several options to obtain independent infor- 
mation on the inclination. For example, we have measured 
vsini ~ 113 km s _1 , which is about 3 times smaller than 
the value reported bv iHutchings fc Cramptoiil Jl98lft . We 
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note that our instrumen tal resolution is higher than that of 
iHutchings fe Cramptonl (83 km s _1 versus ~130 km s _1 ), 
which would favour our determination. This is relevant for 
the orbital inclination (assuming that the Be star's spin axis 
and the orbital axis are aligned) because current observa- 
tional constraints indicate that Be sta rs rotate at ~ 0.7-0.8 
times 2 the critical rotational velocity jPorterl fl996l) . In the 
case of the optical component of LS I +61 303 this would be 
^crit -500-550 km s" 1 . Therefore, if the primary is a normal 
Be star, we are viewing the binary at a low inclination angle 
of —15-20° but we cannot rule out that the Be star may 
be rotating at substantially lower speed than « cr it . However, 
the lack of clear shell lines in our spectra, such as Fell A4584, 
points to a rough upper limit i < 60°. A lower limit to the 
inclination of i ~ 10-15° is obtained from the restriction that 
the rotational velocity of the optical companion should not 
exceed 0.9 times the critical value. Perhaps future accurate 
modelling of the spectral line profiles could provide a better 
constraint on the viewing angle of the Be star and thus on 
the orbital inclination. 

Additional information on the inclination might come 
from the geometry of the relatiyistic j et as derived from 
the VLBI images of iMassi et"ai] J2004I) by assuming per- 
pendicularity between the ejection direction and the orbital 
plane. However, because of the possible precession reported 
by these authors, values in the full range i =0-90° could be 
feasible. 

To conclude with the discussion of the orbital inclina- 
tion, it seems that a crude value of ~30° (with an uncer- 
tainty of some 20° ) is the best compromise given the avail- 
able constraints at this point. It is clear that further ob- 
servations will be required to better define this important 
parameter of the system. With the current knowledge, both 
a neutron star and black hole are equally plausible as the 
compact object in the LS I +61 303 system. 

But there is one very relevant parameter that is in- 
deed well-constrained by our orbital solution, and that is 
the phase of the periastron passage. We obtain a value of 
(/'p eri = 0-23 ± 0.02, in good agree ment with that obtained 
bv IHutchings fe Cramptonl dl98lT) (<jf> pe ri = 0.2), when th e 
phase origin is set by the radio ephemeris of lGregorvl |2002). 
A key question is how this periastron passage phase relates 
to the emission variability of LS I +61 303. There have been 
a number of studies of the occurrences of radio outbusts 
in LS I +61 303 and there is general consensus in that the 
radio emission peaks are dis tributed over a wide phase inter- 
val of about 0.45-0.95 (e.g., iParedes. Estalella fc RiuJl99ot 
lOregorv fc Neisbll2002l) . In contrast, X-ray maxima seem to 
occur at somewhat earlier phases of 0.43-0.47, although this 
is based on only two observatio ns. From three well-studied 
radio outbursts, iGregorvl i2002f) was able to determine the 
phases at which the onset of the outbursts occurred and ob- 
tained values in the interval 0.33-0.40. With our new deter- 
mination of the periastron phase, it is now clear that a delay 
exists beween the closest approach of the system's compo- 
nents and the radio and X-ray outbursts. The onset of the 
outbursts appears to take place about 2.5-4 days after peri- 



2 With a likelihood that this might be underestimated due to 
equatorial gravity darkening effects (see Townsend, Owocki &; 
Howarth 2004, Porter & Rivinius 2003). 



astron passage. The delay is likely due to inverse Compton 
energy losses of relativi stic electrons in the radio jets in the 
periastron vicinity fsee lMassll2004bl for details). This now 
well-established result has strong implications on the models 
devised to explain the emission behaviour of LS I +61 303. 
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